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Abstract. The infrared behaviour of the lattice Landau gauge gluon propagator is discussed, combining results from 
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implications for D(0) discussed. 
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INTRODUCTION AND MOTIVATION 

The link between the deep infrared behaviour of the 
gluon and ghost propagators and confinement, has mo- 
tivated a great effort on computing these quantities on 
the lattice. Besides checking gluon confinement criteria, 
another important goal is to compare recent solutions 
of the Dyson-Schwinger equations with lattice results. 
In particular, the scaling solution [1] predicts a vanish- 
ing gluon propagator and a divergent ghost propagator 
at zero momentum. This solution complies with Gribov- 
Zwanziger [2] and Kugo-Ojima [3] confinement criteria. 
On the other hand, the decoupling solution [4] claims that 
a finite and non-vanishing zero momentum gluon propa- 
gator and a tree level like ghost propagator. The value of 
the zero momentum gluon propagator is connected with 
a dynamical generated gluon mass. 

In this paper we report on our current results for the 
Cucchieri-Mendes bounds in SU(3) lattice gauge theory. 

CUCCHIERI-MENDES BOUNDS 

The Cucchieri-Mendes bounds [5] provide upper and 
lower bounds for the zero momentum gluon propagator 
of lattice Yang-Mills theories in terms of the average 
value of the gluon field. In particular, they relate the 
gluon propagator at zero momentum D(0) with 

"<°> = ^IK<°>I. c» 

where d is the number of space-time dimensions, and N c 
the number of colors. In the above equation, A"^ (0) is the 
a color component of the gluon field at zero momentum, 
defined by 



where A"^ (x) is the a color component of the gluon field 
in the real space. D(Q) is related with M(0) by 

(M(0)) 2 < ffl < N d (N? - 1) <M(0) 2 ) . (3) 

In the last equation ( ) means Monte Carlo average over 
gauge configurations. For convenience we will use the 
definition N cc i = Nj(N^ — 1). The bounds in equation 
(3) are a direct result of the Monte Carlo approach. The 
interest on these bounds comes from allowing a scaling 
analysis which can help understanding the finite volume 
behaviour of D(Q): assuming that each of the terms in 
inequality (3) scales with the volume according to A/V a , 
the simplest possibility and the one considered in [5], an 
a > 1 for (M(0) 2 } clearly indicates that D(0) -» as the 
infinite volume is approached. In this sense, this scaling 
analysis allows to investigate the behaviour of D(Q) in 
the infinite volume limit. 

For the SU(2) Yang-Mills theory [5], the results show 
a D(Q) — for the two dimensional theory, but a D(Q) ^ 
for three and four dimensional formulations. 

RESULTS FOR SU(3) GAUGE THEORY 

We have studied the Cucchieri-Mendes bounds within 
SU(3) lattice gauge theory for three values of the gauge 
coupling: /3 = 6.0 [6, 7], j3 = 5.7 [7], and j3 = 6.2. 

Scaling analysis for /3 = 6.0 

In table 1 we present the lattice setup for j3 = 6.0, 
pointing out the differences to [6, 7]. 

Figure 1 shows the results for the bounds, together 
with the fits to 0}/V a . Assuming this simple scaling 
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FIGURE 1. Cucchieri-Mendes bounds for /3 = 6.0. 

TABLE 1. Lattice setup for p = 6.0. The lattice spacing 
isa = 0.1016(25)fm. 

L 4 16 4 20 4 24 4 28 4 32 4 48 4 64 4 80 4 
L(fm) 1.63 2.03 2.44 2.84 3.25 4.88 6.50 8.13 
#conf. 52 72 60* 56 126 104 120 50 1 ' 



* new ensemble 
1 new statistics 



behaviour, our results for the exponent a support D(0) = 
- see table 2. However, when one assumes a scaling 
behaviour like C/V + (DV~ a , the results support D(0) ^ 
- see table 3. In this sense, a finite and non-vanishing 
value for D(Q) in the infinite volume is not excluded. 

Concerning the fits to co/V a , the reasons for the dif- 
ferences in the values of a reported here and in [5] - and 
therefore on the behaviour of D(Q) in the infinite vol- 
ume limit - are not clear. The simulations use different 
gauge groups. Although there it is generally believed that 
the SU(2) and SU(3) propagators are equivalent for mo- 
menta above 1 GeV [8, 9], a recent direct comparison for 
smaller momenta has shown a measurable difference in 
the infrared region [10]. 

Moreover, the physical volumes used in [5] are much 
larger - up to (27fm) 4 - than the ones used here - up 
to (8fm) 4 . However, the reader should be aware that in 
the SU(2) case the lattice spacing used is about twice the 
lattice spacing considered here. 

TABLE 2. Fits to co/V a using lattice data at /3 = 
6.0. 

to a Xy 

(M(0)) 9.53(36) 0.5255(26) 0.80 
O(0)/V 149 ±10 1.0542(49) 0.63 
N c d{M{0) 2 ) 2927±221 1.0504(54) 0.83 



TABLE 3. Fits to C/V + (0V~ a using lattice data at /3 = 
6.0. 
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LATTICE SPACING EFFECTS IN THE 
GLUON PROPAGATOR 

In order to disentangle possible lattice effects due to 
the use of a different lattice spacing, we carried out 
simulations at j3 = 5.7 and j3 — 6.2. The lattice setup is 
shown in tables 4 and 5 respectively. 

TABLE 4. Lattice setup for j8 = 5.7. The lattice 
spacing is a = 0. 1838(1 l)fm. 

L 4 8 4 10 4 14 4 18 4 26 4 36 4 44 4 
L(fm) 1.47 1.84 2.57 3.31 4.78 6.62 8.09 
#conf. 56 149 149 149 132 100 55* 



TABLE 5. Lattice setup for j5 = 6.2. The lat- 
tice spacing is a = 0.07261 (85)fm. 

L 4 24 4 32 4 48 4 64 4 80 4 
L(fm) 1.74 2.32 3.49 4.65 5.81 
#conf. 51 56 87 99 15 



Some differences have been seen between the gluon 
propagator computed at different lattice spacings for sim- 
ilar physical volumes. An example can be seen in figures 
2 and 3, where the infrared j3 = 6.2 data does not agree 
with data from /3 = 5.7 and 6.0 simulations. These dif- 
ferences deserve further investigations to clarify any pos- 
sible effects due to finite lattice spacing. 



Scaling analysis for /3 = 5.7 and /3 = 6.2 

In what concerns the fits to (o/V a , the analysis of the 
data coming from both sets still supports a vanishing 
D(Q) in the infinite volume limit - see tables 6 and 7. 

Similarly to the case studied before, the lattice data is 
also well described by the functional form C/V + a>V~ a 
- see tables 8 and 9. Although the P — 5.1 case supports 
D(0) ^ 0, for j3 = 6.2 the statistical errors do not allow 
to take any conclusion. In fact, although C = within 
statistical errors, we also get a = 1. For this case, it is 
worth an increase of statistics. 



Renormalized Gluon Propagator - u = 3 GeV 
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FIGURE 2. Comparing the gluon propagator computed 
using different lattice spacings at the same physical volume 
V ~ (4.8/m) 4 . 



FIGURE 3. Comparing the gluon propagator computed 
using different lattice spacings at the same physical volume 
V ~ (6.5/ra) 4 . 



TABLE 6. Fits to C0V~ a using lattice data at TABLE 8. Fits to C/V + coV~ a using lattice data at j8 = 

j8 = 5.7. In order to keep % 2 < 2 > tne 2 6 4 lattice 5.7. In order to keep % 2 < 2, the 26 4 lattice data has been 

data has been excluded. excluded. 
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CONCLUSIONS 

We have studied the scaling behaviour of Cucchieri- 
Mendes bounds using ensembles generated at several 
lattice spacings. Fits of the data to a pure power law in the 
volume strongly support D(0) = 0, but the use of other 
ansatze do not allow to take definitive conclusions. 
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TABLE 7. Fits to ft)V- a using lattice data at TABLE 9. Fits to C/V + coV a using lattice data at /3 = 6.2. 

p = 6.2. Data for M(0) does not include 48 4 . Data for M(0) does not include 48 4 . 
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